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Nonlinear relations among frequencies and phases in modulational instability of circularly polar- 
ized Alfven waves are discussed, within the context of one dimensional, dissipation-less, unforced 
fluid system. We show that generation of phase coherence is a natural consequence of the modula- 
tional instability of Alfven waves. Furthermore, we quantitatively evaluate intensity of wave-wave 
interaction by using bi-coherence, and also by computing energy flow among wave modes, and 
demonstrate that the energy flow is directly related to the phase coherence generation. We first 
discuss the modulational instability within the derivative nonlinear Schrodinger (DNLS) equation, 
which is a subset of the Hall-MHD system including the right- and left-hand polarized, nearly de- 
generate quasi-parallel Alfven waves. The dominant nonlinear process within this model is the four 
wave interaction, in which a quartet of waves in resonance can exchange energy. By numerically 
time integrating the DNLS equation with periodic boundary conditions, and by evaluating relative 
phase among the quartet of waves, we show that the phase coherence is generated when the waves 
exchange energy among the quartet of waves. As a result, coherent structures (solitons) appear in 
the real space, while in the phase space of the wave frequency and the wave number, the wave power 
is seen to be distributed around a straight line. The slope of the line corresponds to the propagation 
speed of the coherent structures. Numerical time integration of the Hall-MHD system with periodic 
boundary conditions reveals that, wave power of transverse modes and that of longitudinal modes 
are aligned with a single straight line in the dispersion relation phase space, suggesting that effi- 
cient exchange of energy among transverse and longitudinal wave modes is realized in the Hall-MHD. 
Generation of the longitudinal wave modes violates the assumptions employed in deriving the DNLS 
such as the quasi-static approximation, and thus long time evolution of the Alfven modulational 
instability in the DNLS and in the Hall-MHD models differs significantly, even though the initial 
plasma and parent wave parameters are chosen in such a way that the modulational instability is 
the most dominant instability among various parametric instabilities. One of the most important 
features which only appears in the Hall-MHD model is the generation of sound waves driven by 
ponderomotive density fluctuations. We discuss relationship between the dispersion relation, energy 
exchange among wave modes, and coherence of phases in the waveforms in the real space. Some 
relevant future issues are discussed as well. 
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I. INTRODUCTION 

In various areas in space and astrophysical environ- 
ment, for instance in the solar wind and in foreshock re- 
gion of planetary bowshocks, parametric instabilities of 
magnetohydrodynamic (MHD) waves are thought to play 
essential roles in generating MHD turbulence. Spacecraft 
observations suggest that, within the MHD turbulence in 
the solar wind and in the earth's foreshock, a large num- 
ber of localized structures are often embedded (Mann et. 
ai, 1994; Dudok de Wit et. ai, 1999; Lucek et. at, 
2004; Tsurutani et. ai, 2005). Also, recent surrogate 
data analysis using magnetic field data obtained by Geo- 
tail spacecraft (Hada et. ai, 2003; Koga and Hada, 2003) 
has revealed that large amplitude MHD waves observed 
in the earth's foreshock are not completely phase ran- 
dom, but are almost always phase correlated to a certain 
degree. Furthermore, the larger the MHD wave ampli- 
tude, the stronger the wave phase correlation, implying 
that the detected phase coherence is a consequence of 
nonlinear interaction among the MHD waves. Since the 
finite phase coherence in the Fourier space corresponds 
to the presence of localized structures in the real space, 
the large amplitude MHD turbulence in the foreshock 



should be regarded as a superposition of random phase 
MHD turbulence plus a large number of localized struc- 
tures. 

Since the MHD waves are dispersive in general, it is 
natural to infer that the localized structures once pro- 
duced would disperse away as time elapses. On the other 
hand, the nonlinearity of the MHD set of equations acts 
as a source for production of the localized structures. 
In fact, inspection of a simple model representing inter- 
action of parallel Alfven wave modes suggests that, the 
phase coherence is generated whenever there is an ex- 
change of wave energy among wave modes in resonance 
(Nariyuki and Hada, 2005). Namely, the localized struc- 
tures are unstable and fade away, but are also contin- 
uously born due to intrinsic nonlinearity of the MHD 
system. Such behavior of the localized structures is typ- 
ically seen in nonlinear evolution of Alfven wave para- 
metric instabilities, in particular, in nonlinear evolution 
of the modulational instability. In a similar context, No- 
cera and Buti (1996) discussed formation of localized 
pulses in a non-dissipative DNLS system, and Hasegawa 
et. al. (1981) observed emergence of organized structures 
in the Korteweg - de Vries equation. The four-wave in- 
teraction scheme was used to explain the emergence of 
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organized structures in the unforced, dissipative DNLS 
(Krishan and Nocera, 2003). In this paper, we discuss 
implications of the presence of finite phase correlation, 
which corresponds to the presence of localized structures, 
generated by the wave- wave interactions in modulational 
instability of circularly polarized Alfven waves, within 
the context of one dimensional, dissipation-less, unforced 
fluid system. We will emphasize the importance of non- 
linear relations among frequencies and phases in identify- 
ing relevant physical process involved in the modulational 
instability. 

The plan of the paper is as follows: In Section 2, we 
review and compare the parametric instabilities in the 
Hall-MHD system and in the DNLS equation. In Sec- 
tion 3, we discuss the relationship between generation 
of solitary waves and nonlinear wave-wave interaction in 
detail, by examining numerically produced Alfvenic tur- 
bulence using the DNLS equation. Our main discussions 
on the phase coherence is presented in this section. In 
Section 4, we discuss phase and frequency features in 
the Hall-MHD system. From the relation between the 
plasma density and the envelope of the magnetic field, 
we examine validity of the so-called quasi-static approx- 
imation used in the DNLS. When the longitudinal and 
transverse wave modes have similar phase velocities, they 
can couple strongly, and as a result the phase coherence 
can be generated. We summarize the results and discuss 
some of the fundamental issues presented in the paper in 
Section 5. 



II. BASIC EQUATIONS AND ANALYTICAL 
STUDY OF PARAMETRIC INSTABILITIES 

It has long been known that circularly polarized ('par- 
ent') Alfven waves are parametrically unstable to gen- 
eration of plasma density fluctuations and 'daughter' 
Alfven waves with the same polarization as the 'par- 
ent'. As early as in the 1960's, it was shown by Galeev 
and Oraevskii (1963), and Sagdeev and Galeev (1969) 
that the circularly polarized Alfven waves are subject 
to parametric decay instability (see Appendix). Later, 
Goldstein (1978) and Derby (1978) derived the disper- 
sion relation of the decay instability of the circularly po- 
larized Alfven waves using ideal MHD equation set. Roles 
of the dispersion effect was investigated using the Hall- 
MHD (two-fluid) set of equations by Wong and Goldstein 
(1986), Longtin and Sonnerup (1986), and Terasawa et. 
al. (1986). The Hall-MHD equations are 
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where the density p is normalized to the initial uniform 
density poj the magnetic field b to the background con- 
stant field magnitude box, the velocity v to the Alfven 
velocity defined by po and b 0x , and the pressure p to the 
ambient magnetic pressure. Time and space are respec- 
tively normalized to the reciprocal of the ion cyclotron 
frequency and the ion inertial length defined using the 
background quantities. 

In this paper, all the physical quantities are assumed 
to be dependent only on one spatial coordinate (x). Then 
the governing equations become 
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where b — b y + ib z and v = v y + iv z are the complex 
transverse magnetic field and velocity, respectively, and 
u = v x is the longitudinal velocity. For simplicity, we 
assume the equation of state to be isothermal, p = (3p, 
where a constant (3 is the squared normalized sound wave 
speed. 

In the following, linear perturbation analysis is used 
to discuss parametric instabilities of parallel propagating 
Alfven waves, using the DNLS equation (Section 2.1) and 
the Hall-MHD equations (Section 2.2). 



A. Linear perturbation analysis of the DNLS 
system 

By applying a quasi-static approximation in which hy- 
drodynamic nonlinearities and steeping are weak (i.e., 
variation of the plasma density is caused only by mag- 
netic ponderomotive fluctuations), Rogistcr (1971) de- 
rived a kinetic equation describing the long time evolu- 
tion of the Alfven waves. Essentially the same equation 
was subsequently obtained starting from the two-fluid set 
of equations (Mj0lhus, 1974; Mio et. al, 1976; Spangler 
and Sheerin, 1982; Sakai and Sonnerup, 1983; Mj0lhus 
and Hada, 1997). The fluid version, now known as the 
DNLS, reads 
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n\ where normalizations are the same as in (l)-(4), a/C^ = 
Cf /4(Cf - Cf), 2(i/C'a = c/bj pl is the ion inertia disper 



sion length, and Ca, Ci, C s are the Alfven, intermediate, 
and sound speeds, respectively. In the DNLS equation, 
b = —v, and 
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where = C% / C\ (the quasi-static approximation) . 

The DNLS describes evolution of weakly nonlinear, 
quasi-parallel propagating, both right- and left-hand po- 
larized Alfven waves (or " magnetosonic" and "shear 
Alfven" waves), which are nearly degenerate. Modula- 
tional instability is driven unstable for the left-hand po- 
larized waves when (3 < 1, and for the right-hand po- 
larized waves when (5 > 1 (Mj0lhus, 1976; Spanglcr and 
Sheerin, 1982; Sakai and Sonnerup, 1983), although pres- 
ence of resonant ions significantly alter the above condi- 
tions, especially for [3 > 1 and moderate to large ion 
to electron temperature ratio (Rogister, 1971; Mj0lhus 
and Wyller, 1988; Spangler, 1990; Medvedev et. at, 
1997). The DNLS is known to be intcgrable under various 
boundary conditions (Kaup and Newell, 1978; Kawata 
and Inoue, 1978; Kawata et. ai, 1980; Chen and Lam, 
2004). In this paper, we restrict our discussion to the 
case (3 -C 1, since only within this regime the use of the 
fluid version of the DNLS is justified, unless the ion to 
electron temperature ratio is extremely low. 

By re-scaling of the variables, t — > [it jo? and x — > 
p,x/a, (9) is reduced as 
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A parallel, circularly polarized, finite amplitude Alfven 
wave, bp — bo exp(i(u>ot — k^x)), where &o, ^o, and fco 
are real, satisfies (11) with the dispersion relation for the 
parent wave, u>o = b^ko + ^o- Now we superpose small 
fluctuations to b p and write 



b = b p + e M b n exp(i(u) n t-k n x)), (12) 

n— — oo 

where k n = k + nK, ui n = uj + nQ, and n is an integer 
(except for 0). Here, SI and K are the frequency and 
the wave number of longitudinal perturbations (quasi- 
modes), u = p. At order of e ' =*= 1 ' , we obtain linearized 
equations 

Q = 2K(b 2 a + ho) ±K^K 2 + 2b 2 a k a + bl (13) 



If < \K\ < boy/— 2ko — 6q, the system exhibits the 
modulational instability. The wave number correspond- 
ing to the maximum growth rate is K max = b (—k — 
&Q/2) 1 / 2 . When the system is modulationally unstable, 
the frequency of sideband modes is obtained from the real 
part of (13) as uj±i = 2v c f,ok±i — u>q, where v^q — cog/kg. 
Therefore, the dispersion relation of the sideband modes 
appears as a straight line, since the dispersion term is 
cancelled by the nonlinear term. In other words, the 
nonlinear effect and the dispersion effect are balanced. 
We will encounter again such a dispersion relation in nu- 
merical analysis in Section 3 and 4, where more detailed 
discussions will be given. 



B. Linear perturbation analysis of the Hall-MHD 
equation 

Finite amplitude, circularly polarized Alfven wave in 
the form below is an exact solution to (5)- (8), 



b p = b cxp(i(uj t - k x)), 
v p = v exp(i(u> t - k x)), 



(14) 
(15) 



where v = -bo/v^o, v^q = uj /k , p = 1, and u = 0, 
together with the dispersion relation 



u 2 = fcg(l + w ). 
Now we add small fluctuations and write 
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where k n — k + nK, uj n = uj + nO, n is an integer, and 
c.c. represents the complex conjugate. At order of e^ 1 ', 
we obtain the following linearized equations 
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where we write the subscripts (±1) as (±) for brevity. 
Combining these equations above, we obtain, 

(n 2 -pK 2 )L + L_ = ^°(L + R_+L_R + ), (27) 



where 



L±=u>±- fe ± (l + u>±), 



(28) 



The dispersion relation (27) has been obtained and 
solved numerically (Wong and Goldstein, 1986; Longtin 
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FIG. 1: Dispersion relation of parametric instabilities in the 
Hall-MHD system. Shown are the real (red lines, left scale) 
and imaginary (blue lines, right scale) normalized frequencies 
plotted versus the normalized wave number, K/ko. Parame- 
ters used are, fco = —0.5, /3 = 0.5, and bo = 0.4. For com- 
parison, the dispersion relation of sideband wave modes and 
the sound wave (zero parent wave amplitude) is superposed 
as black solid and dashed lines, respectively. Forward and 
backward lower (upper) sideband waves are labeled / — (+) 
and b — (+), and the sound wave is labeled s. 



and Sonnerup, 1986; Terasawa et. ai, 1986; Vinas and 
Goldstein, 1991; Hollweg, 1994; Champeaux et. ai, 
1999). Fig. 1 shows the real and imaginary frequencies 
of the longitudinal waves when f} = 0.5, bo = 0.4, and 
fco = 0.5. Under this particular set of parameters, the 
modulational instability (0 < K/ko < 1) has a growth 
rate larger than (or almost comparable with) the beat 
instability around fc = 1.3. 

It is instructive to look at (23) and (24) in order to 
discuss how the side band mode frequencies are deter- 
mined. The r.h.s. of these equations represent three 
basic processes which determine the transverse magnetic 
field; the linear response (k±v±* > ), the v x b nonlinear- 

ity (k±b ui/2), and the Hall effect — b k k±p/2). 

Hereafter we will call these terms as "linear", "nonlin- 
ear", and the "dispersive" terms, respectively. For a 
given set of parent wave parameters (&o and fco), (3, and 
K, one can solve (27) to find £1 associated with the 
modulational instability (eigenvalue), together with ra- 
tios among variables, p\, u±, b± , and v± (eigenvector). 
By varying K in both positive and negative regimes, the 
dispersion relation can be visualized by plotting the real 
part of lo = ujQ+n as a function of fc = ko+K, as shown in 
Fig. 2 as a solid red line. Superposed in the figure are the 
plots of Lo Un = (k±v±^)/b±\ u>„i = k±b Ui/2b±\ and 

^disp = fc± — b k k±p/2b±^ versus fc, which are frequency 
contribution from the three effects introduced above. By 
definition, lo = uo lin + to n i + Lo disp . 

From the figure we notice that both w and ioun al- 
most linearly depend on fc, and thus sum of Lo n i and LOdisp 
should also depend linearly on fc, suggesting that the non- 
linear effect and the dispersion effect are balanced, as 
they were for the DNLS equation discussed previously. 
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-I. G t ^— , . 

FIG. 2: The three frequency contributions, uonn, u) n i, an d 
^disp, plotted versus fc, using the same set of parameters used 
for Fig. 1. See text for explanation of how to make the plot. 
Both Re[tj] (red solid line) and uiu n are linearly dependent on 
k, and so should be the sum of ojdisp and co n i, suggesting that 
the dispersion and the nonlinear effects are cancelled each 
other. 



In other words, the mismatch of the resonance condition 
due to the dispersion effect is reduced by nonlinear inter- 
action among the Fourier modes, just like in the case of 
the DNLS (Section 2.1). 



III. PHASE COHERENCE IN THE DNLS 
EQUATION 

In this Section, we discuss the relationship between 
generation of solitary waves and the wave-wave inter- 
action. In Fig. 3 we show numerical time integration 
of (11) under periodic boundary conditions, where the 
envelope (|6|) is plotted in the phase space of time (t) 
and space (x). As initial conditions, finite amplitude, 
left-hand polarized, monochromatic Alfven waves with 
the wave amplitude bo = 0.4 and the wave mode num- 
ber to = mo = —11 are given, superposed with a very 
small amplitude white noise with < |6 no ise| 2 > 1 / 2 = 10~ a 
within the range of —256 < to < 256, where the bracket 
denotes spatial average over the simulation system. In 
the above, the wave number fc is related to the mode num- 
ber to as fc = 27TTO./L, where the system size is L = 256. 
We have used the convention that to and fc positive (neg- 
ative) represents the right- (left-) hand polarized waves. 
For the numerical computation, we have employed the ra- 
tionalized Runge-Kutta scheme for time integration and 
the spectral method for evaluating spatial derivatives. 
Number of grid points used for this run is 2048. 

Since in Fig. 3 we are plotting the magnetic field en- 
velope, which is constant for the monochromatic, circu- 
larly polarized waves, at the beginning not much wave 
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FIG. 3: Time evolution of envelope |6| in the DNLS model 
(11) with periodic boundary conditions. Initial conditions are 
given as a superposition of finite amplitude, left-hand polar- 
ized monochromatic waves, and very small amplitude white 
noise. 



FIG. 4: Time evolution of the power spectrum (in logarithmic 
scale), log \ b m \, plotted in the phase space of the wave mode 
number (m) and time. The parent mode is given initially 
at mo = —11. Positive (negative) m corresponds to right- 
(left-) hand polarized waves. The wave number is given by 
k — 2irm/L, where L = 256 is the system size. 



activities are apparent. However, starting from t ~ 200, 
modulation of the envelope becomes increasingly more 
evident, and around t ~ 300 a series of solitary waves is 
created. The number of solitary waves is decided by the 
wave number that has the largest growth rate in (13). 
After t ~ 300 we observe a complex behavior of solitary 
waves: they appear, propagate, disappear, and interact 
with each other. 

Corresponding time evolution of the power spectrum is 
shown in Fig. 4. During < t < 200, the parent wave en- 
ergy (m = —11) is gradually transferred to the side-band 
daughter waves (mainly, m = —4 and m = —18) through 
the modulational instability. Later on, increasingly more 
waves at different mode numbers are generated due to 
coupling among finite amplitude waves, and also due to 
the modulational instability of the daughter waves, as 
can be seen as widening of the power spectrum in the 
m-space. Around t ~ 300, the width of the spectrum 
is maximized, corresponding to the appearance of the 
solitary waves. When the solitary waves disappear, the 
power spectrum becomes narrow again, representing the 
uncertainty principle, i.e., the width of the wave packet 
in the real and Fourier space are inversely proportional 
to each other. We note that solitary waves disappear al- 
most completely at certain time intervals like t ~ 900. 
This is reminiscent of the presence of (near) recursion of 
the DNLS equation with periodic boundary conditions. 

Next we show that distribution of wave phases and 
frequencies are closely related to the generation and dis- 
appearance of the solitary waves. We write the Fourier 
transformed magnetic field = \bu\ expi<j>k, and discuss 
the relation between behavior of solitary waves and cor- 
relation among wave phases, 4>k- A method to evaluate 
the wave phase coherence has been proposed (Hada et. 
al, 2003; Koga and Hada, 2003), which we briefly ex- 
plain below: Suppose we have a data containing waves, 
Borg(x), where 'ORG' stands for 'original data'. In the 
present analysis, this is a snapshot of simulation data 
evaluated at certain fixed time. From Bqrg{x) we make 



a phase randomized surrogate (PRS), Brrs(x), and a 
phase correlated surrogate (PCS), Bpcs( x ), by shuffling 
and making equal all the phases, respectively. Then we 
compute the " phase coherence index" , 

Cf, = (Lprs — Lorg)/ (Lprs - Lpcs), (30) 

where 

L* = J2\ B *( X + S )- B *(x)\ (31) 

X 

is the first order structure function evaluated for data (*), 
and the asterisk stands for either ORG, PRS, or PCS. In 
the above, S is an external parameter representing the 
coarsing scale, which we choose to be the grid size for 
the present analysis. When is close to 0, wave phases 
are almost random, while being close to unity sug- 
gests that phases are almost completely coherent. Fig. 
5 shows time evolution of for the run shown in Fig. 
3. Apparently, the appearance and disappearance of soli- 
tary wave trains correspond to the increase and decrease 
of . This is easily understood since the solitary wave- 
form may be produced by equal-phase superposition of 
many waves with different wavelengths. 

Now we discuss how phase coherence is generated by 
the modulational instability in the DNLS. Since the non- 
linearity of the DNLS equation is cubic, it follows that, 
the basic nonlinear interaction in this model is the four 
wave resonance. (On the other hand, if we regard \b\ 2 
as a quasi-mode as represented in the static approxima- 
tion (10), the interaction can be viewed as the three wave 
resonance among two Alfven waves and the quasi-mode.) 

We have the resonance relation among the wave num- 
bers, 

fci + fc 2 = k 3 + h, (32) 

and a similar relation should hold for the wave frequen- 
cies as well. While the resonance condition of the wave 
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FIG. 5: Time evolution of the phase coherence index, Cp. 
When C4, ~ 0, the wave phases are almost random, while 
C$ ~ 1 represents that the wave phases are almost completely 
correlated. 



frequencies usually has the mismatch due to the disper- 
sion effect in the DNLS, as we saw in Section 2, the dis- 
persion effect and the nonlinear effect cancels each other, 
and so there is no mismatch of frequencies when the wave 
modes are coupled. 

Let us define the relative phase among the four waves 

as 



0(k) = 4>ki + 4>k2 ~ 4>k3 - <t> 



(33) 



where k = {ki, ki, &3, fei} represents the quartets of wave 
modes in resonance. Small temporal change of 0(k) im- 
plies that the phase coherence between the four waves 
is strong (locked), because this temporal change corre- 
sponds to the mismatch of the resonance condition of 
frequencies. 

From (11) and (33), we have 
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where the summation is to be taken over all the combi- 
nations of (k2 7 k 3 ,k 4 ) which satisfy (32). Equation (34) 
indicates that the evolution of wave energy of a certain 
mode is determined by 'energy flow' exchanged between 
the quartets, 



F(k) ^MMMIMsin^k). 



(36) 



The direction of the energy flow is determined by the sign 
of sin#(k) and k%. 

Fig. 6 is a scatter plot in the phase space of -F(k) and 
\dd(k)/dt\ (= |0(k)|), in which a single dot is plotted at 
every time step based on the simulation run shown in 
Fig. 3. The set of wavenumbers, kj — 2-Krrij/L, with 
7713 = TO4 = mo = — 11, mi = —15, rri2 = —7 are chosen 
in such a way that the waves satisfy the resonance con- 
dition, and also that the modes m = m\ and m-2 are the 



two daughter waves driven unstable by the parent wave 
at m = mo. It is seen that exchange of the wave energy 
among the four wave modes is enhanced (reduced) when 
the relative phase is close to constant (varies rapidly) 
in time. The same tendency is seen in any quartets in 
resonance with different choice of the wave numbers. Us- 
ing this result, we can now interpret the time evolution 
shown in Fig. 3 and 4 in detail. Around t ~ 200, due 
to the modulational instability driven by the finite am- 
plitude parent wave, a pair of side band waves (daughter 
waves) appears. At this stage, phase coherence is gener- 
ated only within a small number of resonant quartets, 
and the growing sideband waves are restricted within 
k < (when (3 is low), because of the derivative nonlinear 
term in (11): if a quartet includes some of the right-hand 
polarized wave modes (k > 0), this quartet is "stable" in 
the sense that the exchange of energy among them stays 
only at a fluctuation level. Around t ~ 300, the side- 
band waves are sufficiently large, and quartets of higher 
harmonic sideband wave modes (including k > modes) 
become unstable. This broadening of the power spectrum 
(as seen in Fig. 4) corresponds not only to generation of 
solitary waves (Fig. 3), but also to the generation of the 
phase coherence (Fig. 5), since a large number of quartes 
undergo significant nonlinear interaction among them as 
evidenced in Fig. 6. 

The interpretation above is justified by inspection of 
the dispersion relation directly computed from the sim- 
ulation run. By computing <fik using (35) and plotting it 
versus k\, the dispersion relation in the simulation sys- 
tem can be visualized for any given time. In Fig. 7 we 
show time evolution of the dispersion relation, at t — 0, 
t = 200, and t — 550. Initially, almost all the wave modes 
have extremely small amplitude except for the parent 
wave, and thus they are located around the linear disper- 
sion relation, u> = b^k + k 2 . As many of the wave modes 
acquire finite wave amplitude via the wave-wave interac- 
tion, they tend to be aligned on straight lines (two line 
segments exist in (b): notice a little difference in slopes of 
lines for k positive and negative) . The slope of the line 
corresponds to the propagation speed of a soliton, and 
the interval of k which are aligned on the line segment 
represent the waves the soliton is composed of. At later 
time (c), it is more evident that the wave modes are con- 
centrated on a number of different line segments, with 
different ranges of k, and with different slopes. These 
line segments do not necessarily go through the origin, 
since they represent the four (but not three) wave inter- 
action. Consequently, even though there can be many 
solitary waves (coherent structures) with different prop- 
agation speeds, the resonance condition among the four 
waves is well satisfied, i.e., the frequency mismatch is 
small, and thus a large flow of energy exchange among 
the modes is expected. 

The nonlinear interaction among different waves can 
be inspected through evaluation of the so-called bi- 
coherence (Dudok de Wit et. al., 1999; Diamond et. al., 
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FIG. 6: Relation between F(k) and |0(k)| (see text for de- 
tailed explanation of the plot and the variables). The set of 
wave numbers chosen for the plot are kj = 2nm,j/L, where 
L — 256 is the system size, and m.3 = 7714 = mo = —11, 
mi = —15, 7712 = 7713 + 7714 — mi = — 7. The figure sug- 
gests that the exchange of wave energy between the sites is 
enhanced (reduced) when the relative phase is almost con- 
stant (varies rapidly) in time. The unit of the vertical axis is 
degrees per unit time. 



2000), 



bc(ki,k 2 ,k 3 ) 



< Pkibk2b* k3 > 



\/< \pkib k 2\ 2 >< \bksl 2 >' 



(37) 



Presence of finite wave power at ki, k 2 , k 3 , together with 
be ~ 1 suggests that the waves are in resonance. Fig. 8 
compares time evolution of averaged bi-coherence index, 
60 = TV -1 be, total absolute energy flow among all the 
quartets, F to t — N^ 1 l^(k)|, and the phase coherence 
index, Cs. On computing b^ and F tot , all the combina- 
tions of wave numbers (N different ways) are exhausted. 
Apparently variations of the three quantities agree well 
to each other, suggesting that the phase coherence is gen- 
erated as the wave energy is exchanged among resonant 
quartets. 

Finally, we make a remark on the direction of the en- 
ergy flow. From Fig. 4 we see that the parent wave mode 
(k = ko — 2inno/L) remains to be the dominant one 
throughout the simulation run. Therefore, among the 
many quartets of wave modes in resonance (fci, k 2 , ks, k^) 
with ki+k 2 = k 3 + k4, the most dominant quartet (which 
makes (36) the largest) is the one with k 3 = fc 4 = fcrj. The 
resonance condition of wave number for this quartet is, 



k + k a = ki + k 2 



(38) 



In fact, the set of wave numbers used for Fig. 6 was 
chosen so that (38) is satisfied. In Fig. 9 we have plot- 
ted sin#(k) in the phase space of mi and time, with 
= ki = ko and k 2 determined by (32). The plot shows 
the direction of the energy flow between the parent wave 
and the daughter waves. For example, associated with 
the first formation of solitary wave rows (t ~ 200 — 350) , 
sin#(k) < for k 1 < and sin#(k) > for k\ > 0. 
Namely, in both regimes of k\ > and k\ < 0, the en- 
ergy flow defined in (36) is positive, i.e., the energy is 




FIG. 7: The wave frequency <f>k (red circles) versus k evalu- 
ated at (a) t = 0, (b) t = 350, and (c) t = 550. At t = 0, all 
the wave modes except for the parent wave have very small 
wave amplitude, and thus they are located around the linear 
dispersion relation, lo = b$k + k 2 (blue line). However, as the 
wave modes grow via wave-wave interaction, they tend to be 
aligned along straight lines. The slope of the line corresponds 
to the propagation speed of a soliton, and the interval of k 
which are aligned on the line segment represent the waves 
the soliton is composed of. Among the quartet of waves on a 
single line, the frequency mismatch for the resonance condi- 
tion |#(k)| is small, and consequently, energy exchange among 
them F(k) is large, as seen in Fig. 6. 



transferred from the parent wave to the daughter waves. 
During the time interval when solitary waves disappear 
(t ~ 350 — 400), we observe the opposite, i.e., the energy 
flows back from the daughter to the parent waves. 



IV. NONLINEAR RELATION AMONG PHASE 
AND FREQUENCY: HALL-MHD EQUATIONS 

In this Section, we discuss the modulational instability 
in the Hall-MHD equations (Eqs. (5)-(8)), and compare 
the results with that in the DNLS. 
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FIG. 8: Time evolution of C<p (red area), b$ (blue area), and 
F to t/W~ 6 , (black solid line). Variation of Cp well corresponds 
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FIG. 9: Evolution of sin#(k), plotted in the phase space of 
the wave mode number (mi) and time. 

Numerical time integration of the Hall-MHD set of 
equations, (5)-(8), is performed, in a similar way as for 
the DNLS equation discussed in the preceding section. 
As initial conditions, finite amplitude, left-hand polar- 
ized, monochromatic, parallel Alfven waves are given, 
with the wave amplitude 60 = 0.4 and the mode number 
m = mo = —20. Superposed with the parent wave is a 
small amplitude white noise with < |a„ oise | 2 > 1 / 2 = 10~ 5 , 
where a represents any variable in (5)- (8) within the 
range of —128 < m < 128. The system size is L = 256, 
in which there are 1024 grid points. We chose f3 = 0.5 
so that the modulational instability growth rate is the 
largest among all the parametric instabilities possibly 
driven unstable. The parameters specified here are essen- 
tially equivalent to the ones used in the previous section. 

Fig. 10 shows time evolution of |6| plotted in a way 
similar to Fig. 3, and Fig. 11 represent time evolution of 
log \ b m \ and log \p m \, which is the power of the transverse 
and the longitudinal waves, respectively, as a function of 
the wave mode number, m. 

From comparison of the DNLS and the Hall-MHD sim- 
ulation runs, we dicuss validity of various approximations 
assumed in the DNLS. In particular, we look at the fol- 
lowing: (a) the quasi-static approximation, (10), (b) con- 



FIG. 10: Time evolution of envelope |6| due to Eqs.(5-8) 
with periodic boudary conditions. Initial conditions are given 
as a superposition of finite amplitude, left-hand polarized 
monochromatic waves, and very small amplitude white noise. 

servation of the magnetic field energy alone, and (c) the 
assumption of constant plasma density in the dispersion 
term (i.e., the Hall-term, d/dx(p~ 1 db/dx) 1 is simplified 
by assuming p to be constant). These approximations are 
justified as long as the power of the longitudinal modes 
(Fig. 11(b)) is much less than that of the transverse 
modes (Fig. 11(a)). Around t ~ 350, the power spectrum 
of sideband wave modes and the daughter wave modes 
begin to grow. To check the validity of the approxima- 
tion in (10), we evaluate the cross-correlation function, 

c(A,r) = - <l^ + V + r)|^,t)> 

V< |6(x + A,t + r)| 4 >< (p{x,t)) 2 > 

where we simply let A = r = in the present analysis. 
The result is plotted as a solid line in Fig. 12(a). When 
< t <^ 350, relation (10) is well held, but around 
t ~ 350, the cross correlation starts to decrease. Fig. 
12(b) shows time evolution of the magnetic field energy 
(broken line), parallel kinetic energy (solid line), and per- 
pendicular kinetic energy (dotted line) . The decay of the 
correlation function occurs at the same time as the mag- 
netic field and the perpendicular kinetic energy both de- 
crease. Around t ~ 500, inverse cascade begins to take 
place (discussion of the inverse cascade in the DNLS can 
be found in Krishan and Nocera, 2003). 

The validity of the static approximation has been ex- 
amined by several authors. By performing hybrid (ki- 
netic ions + an electron fluid) simulations, Machida et. 
aZ.(1987) argues that the system automatically adjusts 
itself to a state described by the static approximation 
when ion kinetic effects are included. However, as time 
elapses, longitudinal mode grows and the static approxi- 
mation is soon violated. Spangler (1987) concluded that 
the breakdown of the static approximation takes place 
in association with rapid evolution of wave packets by 
pondcromotive force. Combining (5) and (6), we have 

d 2 p a d 2 p d 2 ,\b\ 2 2 , 
This equation describes propagation of sound waves in 
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FIG. 11: Time evolution of power spectrum (in logarithmic 
scale), log \b m \ and log \p m \, plotted in the phase space of the 
wave mode number (m) and time. The parent mode is given 
initially at mo = —20. Positive (negative) m corresponds to 
positive (negative) helicity. 



a presence of a source term (r.h.s.), which consists of 
the ponderomotive force, \b\ 2 /2, and the advection effect, 
pu 2 (Spangler, 1987). The latter represents the "self- 
coupling" of longitudinal waves. 

In order to clarify the roles of the " self-coupling" term, 
which is not discussed in the analysis of Spangler (1987), 
we have run numerical experiment in which the advec- 
tion terms are artificially removed from the Hall-MHD 
set of equations. The broken line in Fig. 12(a) depicts 
time evolution of the cross correlation function evaluated 
for the run without the advection terms, using the same 
simulation parameters as the previous run. The time 
evolution for both runs are essentially the same, suggest- 
ing that the process playing the most dominant role in 
violating the static approximation is the ponderomotive 
force, as suggested by Spangler (1987). The self-coupling 
term cannot be neglected, however, since its magnitude 
becomes comparable to other terms as the longitudinal 
wave modes grow. 
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FIG. 12: Time evolution of (a) the cross-correlation, c(A = 
0, t = 0) evaluated for the run shown in Fig. 10 (solid line), 
and c(0, 0) evaluated for a run in which the advection term is 
artificially removed from the Hall-MHD set of equations (bro- 
ken line), (b) magnetic field energy (broken), perpendicular 
kinetic energy (dotted), and parallel kinetic energy (solid), 
plotted versus time. Around t ~ 350, together with the decay 
of the cross-correlation function, the magnetic field and the 
perpendicular kinetic energy decrease and the parallel kinetic 
energy increases. 



In a similar way described before to plot the dispersion 
relation of nonlinear Alfven waves in the DNLS directly 
from the simulation run, we have computed the <j) p k and 
(j) b k , using (5) and (8), respectively. Corresponding to ap- 
pearance of solitary waves and localized structures, the 
frequency distribution ((/>£ and <jr k ), i.e., the dispersion 
relation, appears along a straight line, as shown in nu- 
merically evaluated dispersion relation, Fig. 13. Since 
these waves propagate in the same direction as the par- 
ent wave, the frequency distribution along a straight line 
is usually satisfied at the origin. Furthermore, Fig. 13 
shows that the longitudinal and transverse wave modes 
have a similar slope, which is the velocity of the struc- 
tures composed of these waves. Thus, the three wave 
resonance condition between the transverse wave modes 
and the longitudinal wave modes is approximately satisi- 
fied. 

Fig. 14 displays time evolution of the wave frequency 
for the parent wave (m = —20, solid line), and some 
of the daughter waves (m — —8, dotted line) (m = —3, 
broken line). When < t <^ 500, the parent wave 
mode is dominant, and its frequency stays almost con- 
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stant. Around t ~ 500, the inverse cascade begins to take 
place, and the parent wave frequency starts to be modi- 
fied. On the other hand, the frequencies of the lower side- 
band wave modes remain almost constant (broken line). 
Thus, the phase velocity of the modes with the maxi- 
mum power is close to the velocity of the structures, and 
it corresponds to the slope of the dispersion relation in 
Fig. 13. 

Now we examine the approximation (c) we listed before 
in the DNLS, i.e., the assumption of constant density in 
the Hall term. In order to do so, let us explicitly write 
down the phase relation of (8), 



^disp 



dt 



■ E 

k=kl + k2 



\uki\\bk2\ 



COB(0£ 



Irtt h 



k^k 



■■•>'■ £ 

fc=fc3+fc4 



\b k3 \\V M \ 



cos 



=kl+k2/ 



qbbV \ 
7 k=k3+k4> 



I Oft 



(41) 



(42) 



(43) 



(44) 



where V = p^ 1 and 9 k ^ l+m — ^ — (f^ — cf)^. In the linear 
perturbation analysis in section 2.2 we have called cu n i, 
LOdisp, and uJun as contribution to the wave frequency 
via 'nonlinear', 'dispersion', and 'linear' effects, respec- 
tively. The dispersion effect is not simply k 2 , but is 
k 2 — bokoK p/2/b± , where the latter term arises due to 
longitudinal perturbation. Just like in the DNLS model, 
the balance between the 'nonlinear' and the 'dispersion' 
terms also exist in the Hall-MHD equations. 

Fig. 15 shows uj n i, uJdisp, and w; m evaluated at t=350, 
500, and 720. We see that uu n is linear in k: this 
may be regarded as self-generation of the (dispersive) 
Walen's relation, which is the relation between the trans- 
verse magnetic and the velocity field satisfied for AlfVenic 
fluctuations and rotational discontinuities {e.g., Landau 
and Lifshitz, 1962). For example, the normalized mag- 
netic and the velocity field, bo and vq in (14-15), satisfy 
t>o = T^o/f^o, where the signs represent parallel/ anti- 
parallel propagation with respect to the background mag- 
netic field, and the phase velocity u^o is determined from 
the dispersion relation (16). The term cos(<fi b k — in 
the r.h.s. of (44) takes the value of either 1 or -1 (so 
that LOii n /k = ±|ffc|/|6fc|) in association with broadening 
of the power spectrum. Also, from Fig. 13 and 15, we 
find numerically that ujn n /k ~ k/uj (corresponding to the 
appearance of near straight lines in these plots). There- 
fore we find that the system approaches automatically 
to a state where k/cu = ±\vk\/\bk\ (Walen's relation) is 
automatically fulfilled. Furthermore, since uu n is almost 
linear in k in Fig. 13, so is the sum of UJdisp and uj n i as 




FIG. 13: Numerically obtained dispersion relation at (a) 
t=350, (b) t=500, (c) t=720. The blue line indicates the lin- 
ear dispersion relation, and the red circles and black crosses 
show the frequency evaluated from the temporal variation of 
the magnetic field mode (</>y and that of the density modes 
respectively. Superposed (green) line has a slope equal 
to the phase velocity of the wave mode which has the max- 
imum power. Corresponding to the appearance of solitary 
waves and localized structures, the wave power tends to be 
aligned along a straight line. 



(41) suggests, i.e., the contribution to the wave frequency 
nonlinear in k is cancelled due to balance between the 
dispersion and the nonlinear effects. 

Fig. 16 shows distribution of 4> k (similar to 0?) versus 
k. When < t < 500, these longitudinal fluctuations are 
found to be distributed along the same straight line in the 
dispersion relation, and should be recognized as the pon- 
deromotive density fluctuations, produced via interaction 
of Alfven waves. In the present case, the fluctuations are 
produced in the course of the modulational instability, 
and the two interacting Alfven waves are on the same 
branch with nearly equal wave numbers and the wave 
frequencies, so that the resultant ponderomotive fluctua- 
tions have similar phase velocity as the interacting Alfven 
waves. The phase velocities of the longitudinal fluctua- 
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FIG. 14: Time evolution of the wave frequency (<fi k ) for the 
waves with m = —20 (solid line), m = —8 (dotted line), and 
m = —3 (broken line). The wave number of these modes are 
k = -0.49, -0.19, and -0.07, respectively. When < t < 500, 
the parent wave mode is dominant, and its frequency is almost 
constant. Around t ~ 500, the inverse cascade begins to take 
place, and the parent wave frequency starts to vary. On the 
other hand, the frequencies of the lower sideband wave modes 
remain almost constant after t > 500. 



tions (p and u in Fig. 16), transverse fluctuations (6 
and v in Fig. 13), and the acoustic wave (solid line in 
Fig. 13), are measured as ~ 0.78, - 0.78, and ~ 0.70, 
respectively. The former two are (almost) equal, and 
are distinct from the last one. Thus we conclude the 
longitudinal fluctuations observed at < t < 500 are 
the ponderomotive density fluctuations, rather than the 
sound wave. Around t ~ 500, however, the dispersion re- 
lations of the longitudinal fluctuaions start to approach 
the sound wave branch. This is a consequence of the 
nonlinear development of the modulational instability, in 
which the ponderomotive density fluctuation acts as a 
driver to excite the sound waves, as represented in (40). 

This excitation of the sound wave is also related to 
particle heating due to the modulational instability of 
Alfven waves. Many authors discussed Landau damp- 
ing associated with envelope modulation of Alfven waves 
(Mj0lhus and Wyller, 1986,1988; Spangler; 1989, 1990; 
Medvedev and Diamond, 1996; Passot and Sulem; 2003). 
However, the ion acoustic wave mode is exclueded, as 
long as the "ponderomotive" expansion is used in scal- 
ing of the system. In the case of the weak nonlinear- 
ity, the "triple-degenerate" expansion makes a similar 
description at j3 ~ 1 (Hada, 1993). The kinetic "triple- 
degenerate" DNLS equation will be useful for elucidating 
the physics of the heating process. The use of numeri- 
cal simulation models (Vasquez, 1995) and/or theoreti- 
cal models (Snyder, 1997; Passot and Sulem, 2004) that 
include the strong nonlinearity and kinetic effects are de- 
sired. 



V. CONCLUSIONS AND DISCUSSIONS 

In this paper, we discussed the nonlinear relation 
among phases and frequencies in the modulational insta- 
bility of parallel propagating Alfven waves in the context 



I- 
1- 

[. 

-I- 


8 w 




- ..*V ' 


■ * 




■r "i K 


-It 


1 




■ 


t -i 


1 — 

it 
1 


1 I 


l~ 
1 

6- 


IW 

■ 




*JP* - - 

V - 


I- 

-J- 




■ 


1 — * — 1 — 

;Jt -1 


D 

Jl 


! i 




FIG. 15: The three basic elements of the wave frequency, 
uju-n (light blue squares), u> n i (blue circles), and tvdisp (green 
crosses), defined in (42)-(44), plotted versus k, evaluated at 
(a) t = 350, (b) t = 500, and (c) t = 720. As discussed in 
Section 2.2, cui in is almost linear in k, suggesting that the dis- 
persive Walen's relation is held. Since the dispersion relation 
(Fig. 13) suggests that the wave frequency, which is a sum 
of the three terms above, is linear in k also, the linearity in 
k of the sum of uidiap and ui n i should be held as well, indicat- 
ing that the dispersion and the nonlinear effects cancel each 
other, and the system is brought to a state where the wave 
resonance condition is easily satisfied. 



of one-dimensional two fluid models, using linear pertur- 
bation analysis and numerical experiments. 

The main results of the present paper may be summa- 
rized as follows: 

(1) Corresponding to generation of solitary waves or 
coherent structures by the wave modulation, dispersion 
relation of the waves tend to be aligned along straight 
segments of lines (Fig. 7 and Fig. 13). The slopes of 
these straight lines correspond to velocities of the solitary 
waves and the structures, and are well estimated as the 
phase (~ the group) velocity of the wave with the largest 
wave power. 

(2) The exchange of wave energy among the resonant 
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FIG. 16: Numerically computed frequency <j>^(ied circles) 
and <j) p k (black crosses) versus the wave number k at (a) t=350, 
(b) t=500, (c) t=720. When < t < 500, the longitudinal 
modes <j) k are distributed around a line slightly (but dis- 
tinctly) above the sound wave branch (solid black line). The 
former may be called the " ponderomotive density fluctuation" 
branch. As later times, however, we find that the longitudinal 
wave frequencies are re-distributed around the sound branch. 



wave modes is enhanced (reduced) when the relative 
phase is close to constant (varies rapidly) in time (Fig. 
6). This is a universal feature in a wide variety of non- 
linear systems: e.g., it can also be seen in a single triplet 
model (Appendix). 

(3) The fact that the dispersion relation is essentially 
given as segments of lines in the dispersion relation sug- 
gests that the mismatch in the resonance conditions 
among the coupled wave modes is reduced automatically. 
The waves can efficiently exchange wave energy and gen- 
erate the phase coherence. 

In the DNLS equation (Section 3), increase and de- 
crease of the phase coherence index well corresponds to 
the appearance and disappearance of solitary wave trains. 
However, in the Hall-MHD equations (Section 4), evolu- 
tion of the phase coherence index turns out to be more 



complex. There are two issues to be discussed: the phys- 
ical processes leading to the variation of the phase coher- 
ence index, and the validity of the index. 

In the DNLS, the energy flow among the wave modes 
generates order not only among the wave frequencies but 
also among wave phases (see Appendix, and also Fig. 9), 
and the associated broadening of the power spectrum al- 
ways corresponds to the global energy flow because of 
the balance between the dispersion effect and the non- 
linear effect. However, in the Hall-MHD, the presence 
of the localized structures does not always correspond 
to the broadening of the power spectrum. In Fig. 17 
we plot time evolution of and b^, based on the Hall- 
MHD simulation run discussed before. The broadening 
of the power spectrum correlates well with b^, but not 
quite well with C^. This probably is due to the fact that 
the structures observed in the simulation are not always 
generated by the local energy flow. 

In case of the modulational instability starting from 
a single parent wave plus small amplitude background 
noise, one can evaluate the relative phase among the res- 
onant wave modes and also the energy flow among them, 
and identify clearly time intervals the phase coherence is 
generated (see Fig. 9, and also early part of Fig. 18). On 
the other hand, after the initial stage of the modulational 
instability, wave modes at low wave numbers produced 
via the inverse cascade become dominant, it is no longer 
possible to identify the time intervals in which the phase 
coherence is produced. 

Regarding the phase coherence index, we have to say 
that its physical meaning is still not completely clear, 
although it was demonstrated that the index is a useful 
measure to identify the presence of wave- wave interaction 
(Koga and Hada, 2003). Since the concept of the phase 
synchronization is universal (e.g., He and Chian, 2005), 
it is important to develop a method that can characterize 
the underlying nonlinearity of the system. This discus- 
sion is translated into the explanation or improvement of 
the phase coherence index. 

The relation between the longitudinal and transverse 
modes has important information on the nature of the 
parametric instability. In the dispersion relation, lon- 
gitudinal and transverse modes were found to be both 
linear, with almost the same slope, corresponding to the 
structures generated by modulational instability. On the 
other hand, the longitudinal fluctuations generated by 
the decay instability distribute on the sound wave line 
with some band width (Terasawa et. ai, 1986). Such 
a cross relation among different variables in MHD tur- 
bulence is very important to understand the turbulent 
dynamics. Henceforth, we have to analyze in detail the 
general rules of the relation among these variables in par- 
allel Alfven turbulence generated by parametric instabil- 
ity. Also, it is essential to include the kinetic effects when 
/3 is moderately large. 

It is important to identify the phase coherence brought 
by nonlinear wave-wave interaction, in order to discuss 
underlying physical processes leading to the generation 
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FIG. 17: Time evolution of C<p (red area) and &</, (broken line). 
In association with the broadening of the power spectrum, 
varies accordingly, but Cj, does not. 
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via strong and correlated impulse force rather than via 
sequence of random forces which last for a long time. 

For example, let us consider pitch angle diffusion due 
to the turbulence consisting of superposition of parallel, 
circularly polarized Alfven waves (slab model). Within 
this model, it is well known that the particles cannot dif- 
fuse across 90 degrees pitch angle, within the framework 
of the quasi-linear diffusion, simply because there are no 
waves which can resonate with particles without parallel 
velocities. However, the mirror force can reflect the parti- 
cles quite easily, if there exist solitary waves with a finite 
amplitude variation of the total magnetic field. Further- 
more, in a presence of many of these solitary waves, in- 
dividual particle trajectory may be given as combination 
of quasi-ballistic motion between the solitary waves and 
trapping by one of the solitary waves. Fermi type accel- 
eration is possible (Kuramitsu and Hada, 2000), and fur- 
thermore, depending on relative dominance between the 
ballistic and trapped trajectories, ensemble of these par- 
ticles may appear as either super- or sub-diffusive (Zim- 
bardo et. at, 2000; Carreras et. ai, 2001). In order 
to properly describe time evolution of the ensemble, one 
has to invoke the fractal diffusion formalism (Metzler and 
Nonnenmacher, 1998; del-Castillo- Negrete et. ai, 2004). 

Acknowledgments 



FIG. 18: Time evolution of the relative phase between the 
parent and the other wave modes, sin(<j^, — <f> kl — <f>x). 



of the coherent structures. This is particularly impor- 
tant for analysis of spacecraft data, since there are many 
localized structures in space plasmas which are produced, 
presumably, by processes other than nonlinear wave- wave 
interactions. Dudok de Wit et. al. (1999) and Soucek et. 
aZ.(2003) discussed the so-called Volterra model, a way 
to identify nonlinear wave-wave interactions using the 
satellite and the simulation data, by utilizing the high or- 
der spectrum analyses and weak turbulence theory. This 
method may shed light on understanding the nature of 
wave-wave interactions when the weak and strong tur- 
bulence processes co-exist. In the future, we hope to 
clarify the process of the phase coherence generation via 
the wave-wave interaction in more detail, and to extend 
the Volterra-type models in such a way that the weak 
turbulence approximation is not necessary. 

Finally, we make a remark on energetic particle trans- 
port in the turbulence including the localized structures. 
Since these structures have broad power spectrum in the 
k space, particles with wide range of energy can resonate 
with them. However, this is not the more important side 
of the wave-particle interaction, from the viewpoint of 
the phase coherence - rather, it is in the fact that large 
amount of electromagnetic field and plasma energy is 
concentrated within the solitary wave (phase correlated 
wave/structure), so that it can influence particle motion 
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Appendix: Three wave resonance: A single triplet 

A triplet is known as the most basic element of non- 
linear interaction that can be extracted from various 
phenomena. Decay instability of Alfven waves in space 
plasma is an example. A finite amplitude Alfven wave 
propagating parallel to the ambient magnetic field decays 
into a backward propagating Alfven wave and a forward 
propagating ion acoustic wave. We consider a set of three 
wave equations (3W) discussed by Sagdeev and Galeev 
(1969), 



Ci = 


-iC^C 3 




-tC?C 3 


c 3 = 


-iCiC 2 



where Cj is the normalized complex amplitude of the 
jth mode. The above set of equations is derived via res- 
onance condition among frequencies of each mode, 

u> 3 = u>i + u 2 - (48) 

By introducing the wave quanta (wave action), Nj = 
\Cj\ = €j/uij, where €j is the wave energy of mode j, we 
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FIG. 19: Time evolution of Nj with the initial condition, 
N3 >> JVijiVb. The broken, dotted, and solid lines indicate 
Ni, N2, and N3, respectively. The system is 'unstable' in that 
the 'quanta' originally stored in mode 3 are re-distributed to 
modes 1 and 2 within an instability time scale (~ 3000). 



have the Manley-Rowe relations, 

N 1 +N 3 = const, 



N!~N 2 



const. 



(49) 



(50) 



Now let us write Cj 



Aj exp (i<j>j), with Aj(- 



and 4>j real, then we obtain a relation about the phase 
difference among the triplets (0 = <p 3 — <p2 — <fii), 



fa-<h-<t>i, 
,A X A 2 A 2 A 3 



A 3 A X 



A?, 



A x 



) cos 9, 



cot9—log(A 1 A 2 A 3 ), 



(51) 
(52) 

(53) 



where (f>j is the (nonlinear) frequency and 9 represents the 
frequency mismatch in the resonance condition. From 
the above, the following conservation law is derived 



A\A 2 A 3 cos# = D = const. 



(54) 



The 3W equations at a glance have six degrees of free- 
dom, but if we write them into the real amplitude and 
phase, we find that the equations actually have only four 
degrees of freedom: the real amplitudes A\, A 2 , A 3 , and 
the phase difference 9. Since there are three invariants 
((49), (50), (54)), the 3W equation set is integrable. All 
initial problem of 3W can be solved, and the solution 
can be written using elliptic functions. Behavior of the 
system may be physically interpreted as either stable or 
unstable depending on initial distribution of wave quanta 
to the three eigenmodes (Fig. 19 and 20). 

Let us define the flow of quanta F, which represents 
the interaction between the eigenmodes, 



F = iVi = N 2 = -N 3 = 2A 3 A 2 A 1 sinfl, 



(55) 



and consider its relationship to temporal change of the 
phase difference 9. The sign of sm9 determines the di- 
rection of quanta flow. 



FIG. 20: Same as Fig. 19 except that the initial condition 
here is N2 » Ni,N$. The broken, dotted, and solid lines 
indicate Ni, N2, and N3, respectively. The system is 'stable' 
in the sense that variations of Ni, N2, and ./V3 are small. 
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FIG. 21: Time evolution of iVi, 
that A 1 A 2 A A cos6 = 0. 



N2, N3, and 6 for the case 



As for the initial conditions there are two distinct pos- 
sibilities. The first is that there are some modes without 
any quanta, and the second is that all the modes have 
some finite number of quanta. In the first case, cos 9 = 
always since D = 0. In Fig. 21, we see that 9 jumps 180 
degrees periodically because the direction of the quanta 
flow remain unchanged per half the oscillation period. 

From now on, we only discuss the case where all the 
modes have some finite number of quanta. Since D = 
const ^ 0, we can rewrite F as 



F = 2D tan 9, 



(56) 



and therefore, the condition — n/2 < 9 < it/2 is required 
in order that F be finite. This condition is evident from 
(54). Eq. (56) suggests that \F\ is large when 9 is close 
to either -tt/2 or n/2, while F = when 9 = (Fig. 
22(a)). 

Fig. 22(b) shows time evolution of <j>j. Both fast and 
slow temporal changes of 9 can be seen. Around the point 
where the number quanta is close to take extremal values, 
the temporal change of 9 becomes rapid. If the system is 
unstable, temporal change of the phase difference is much 
faster when the quanta is at local minimum, than when it 
is at local maximum. Fig. 22(c) shows this relationship 
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more clearly, 
we have 



Using (54) and (56), and fa = —D /Nj, 



-A)( 



N 3 



1 



No 



F 



(57) 



where G = AiA 2 A 3 . This equation exactly shows the 
relationship discussed above. 



FIG. 22: (a) The energy flow F (solid line, left axis) and sin9 
(dotted line, right axis) are plotted versus time. When the 
energy is exchanged between the modes, the relative phase 
stays at almost at a constant level. The same tendency is 
observed in systems with multiple degrees of freedom (c/., 
Fig. 9 and 18). (b) Time evolution of fa. The thin, solid, 
and dotted lines indicate fa, fa, and fa, respectively, (c) 
Relation between 6 and F for a (single) triplet. The exchange 
of quanta between the sites is enhanced when the relative 
phase is approximately constant in time, while the exchange 
of quanta is reduced when the relative phase varies rapidly. 



In summary, if the temporal change of the phase dif- 
ference is small (large), quanta flow between the modes 
is large (small). In other words, the stronger (weaker) 
the modes are correlated, are the interaction between 
the modes strong (weak). This relationship between the 
phase coherence and the interaction among cigcnmodes 
is universal. For example, similar relationship is held in a 
system where many triplets are connected, in the DNLS, 
and in the Hall-MHD models. Also, when cos6> = ini- 
tially, we have a situation similar to the case that Nj = 
for one (or more) of the modes, and thus cos = always. 
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